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Abstract 

A Chandra study of pulsar wind nebula around the young energetic pulsar PSR B1509— 58 is presented. 
The high resolution X-ray image with total exposure time of 190 ks reveals a ring like feature 10" apart from 
the pulsar. This feature is analogous to the inner ring seen in the Crab nebula and thus may correspond 
to a wind termination shock. The shock radius enables us to constrain the wind magnetization, a > 0.01. 
The obtained a is one order of magnitude larger than that of the Crab nebula. In the pulsar vicinity, the 
southern jet appears to extend beyond the wind termination shock, in contrast to the narrow jet of the 
Crab. The revealed morphology of the broad jet is coincident with the recently proposed theoretical model 
in which a magnetic hoop stress diverts and squeezes the post-shock equatorial flow towards the poloidal 
direction generating a jet. 

Key words: ISM: jets and outflows — ISM: individual (G320.4— 01.2) — ISM: supernova remnants 
- stars: pulsars: individual (PSR B1509-58) — X-rays: ISM 
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1. Introduction 

Recent X-ray observations revealed that part of pulsar 
wind nebulae (PWNc) have common structures of "Torus" 
and "Jet", as typified by the crab pulsar (Hester et al. 
2002; Weisskopf et al. 2000), the Vela pulsar (Pavlov et al. 
2003; Helfand et al. 2001), and PSR B1509-58 (Gaensler 
et al. 2002) . However the physical mechanisms which form 
the ubiquitous features, especially the jets, are still un- 
clear, because the origin of the jet in the Crab is too com- 
pact to discriminate its structure even with the current 
X-ray observatories. 

The theoretical approaches toward the axisymmctric 
PWNe have been aimed to figure out the Crab nebula. 
The torus in the Crab had been successfully explained by 
one-dimensional MHD simulations (Kennel and Coroniti 
1984a; Kennel and Coroniti 1984b). A young rotation 
powered pulsar is believed to lose its rotating energy via a 
magnetized particle flow from a pulsar, the so-called pul- 
sar wind. Since the pulsar wind is flowing almost at the 
speed of light, it is stalled by a termination shock just 
around the pulsar. When the wind passes through the 
termination shock, the kinetic energy of the wind is con- 
verted into internal energy. Essentially, the post-shock 
flow is decelerated yielding to the conservation law of par- 
ticle flux nvr 2 — const, where n is the number density 
(that is almost constant), v is the flow velocity and r is 
the radius from the pulsar. At the same time, the frozen- 
in condition, rvB = const, induces the compression of the 
toroidal magnetic field B, in which the internal energy 



is efficiently converted into synchrotron radiation forming 
an X-ray torus. 

In practice the magnetic pressure prevents the decelera- 
tion of the post-shock flow, thus the magnetization of the 
pulsar wind is closely related to the geometric structure 
of the PWN. In this way, several methods to calculate 
the magnetization parameter a have been proposed (Rees 
and Gunn 1974; Kennel and Coroniti 1984a; Gallant and 
Arons 1994). Moreover, recent two-dimensional numerical 
studies suggested that the magnetization seems to con- 
trol the formation of the jet (Lyubarsky 2002; Del Zanna 
et al. 2004; Bucciantini 2006). Therefore it is essential to 
evaluate the magnetization parameters to understand the 
formation process of PWNe. 

So far, the shock front of the pulsar wind has been ob- 
served only in the Crab nebula (Weisskopf et al. 2000) 
because of the small angular sizes of PWNe. Nevertheless 
the Chandra X-ray observatory with sub-arcsec-scale spa- 
tial resolution has the potential to detect the wind ter- 
mination shocks in other PWNe. In order to obtain a 
general understanding of the formation process of PWNe, 
we investigate the fine structures around the pulsar PSR 
B 1509— 58 utilizing the Chandra data. 

PSR B1509— 58 is emitting 150 ms radio, X-ray, 
and Gamma-ray pulsation (Seward and Harnden 1982; 
Manchester et al. 1982; Ulmer et al. 1993). From the spin 
parameters, a characteristic age r c = 1700 yr, a spin-down 
luminosity £ = 1.8 x 10 37 ergs s _1 , and a surface mag- 
netic field Bp = 1.5 x 10 13 G have been obtained, making it 
one of the youngest, the most energetic, and highest field 




Fig. 1. (Left) — Overview image of the PWN accompanied by PSR B1509— 58 with 190 ks exposure taken by Chandra. The image 
was smoothed with a Gaussian of <r = 2.0". The white rectangle indicates the FoV of the right panel. (Right) — Close-up image of 
the pulsar vicinity smoothed with a Gaussian of a = 0.5". The inner panel schematically describes the remarkable structures. 



pulsars known (Seward and Harndcn 1982; Kaspi et al. 
1994). The pulsar is accompanied by a bright and large 
PWN with remarkable jet-like features extending to the 
southeast and northwest (Gaensler et al. 2002). 

The pulsar and its PWN appears to be embedded in a 
30 arcmin radio shell G320.4-01.2 (Caswell et al. 1981). 
An HI observation yielded the distance of d = 5.2 ±1.4 kpc 
from the earth (Gaensler et al. 1999). Adopting this dis- 
tance and standard parameters of the ISM and supernova 
explosions yields an age of 6-20 kyr, which is one order 
of magnitude larger than the pulsar's characteristic age 
(Seward et al. 1983). 

Moreover Yatsu et al. (2006) detected proper motions 
of X-ray clumps in RCW 89 which coincide with the north 
radio shell from the 4.3 yr base-line Chandra observations. 
These results imply that the pulsar and the SNR have the 
same progenitor in a cavity. If this is the case the termina- 
tion shock of the pulsar wind can be easily discriminated 
by Chandra. 

The paper is structured as follows. The Chandra ob- 
servations and their results are shown in §2. The imaging 
analysis and spectral analysis arc described in §3 and §4, 
respectively. Finally, we discuss the obtained results in 
§5. 

2. Observations 

The vicinity of PSR B1509— 58 has been observed sev- 
eral times by the Chandra X-ray Observatory. For this 
study, the four most recent observations which are se- 
quentially performed from 2004 to 2005 arc utilized. The 



exposure time of the each observation is about ~ 50 
ks. The aiming point (ACIS-I) was set on (15 h 13 m 55?60, 
— 59°08'08."9). The observational settings are summarized 
in Table 1. 

We have performed the reduction process for the each 
data set using "acis_event_process" in CIAO version 4.0. 
Although these four observations are performed in VFaint- 
Mode, we applied standard background cleaning using 
3x3 event island, because the VFaint-Mode background 
cleaning process tends to misjudge a certain amount of 
source photons to be background events. 

Figure 1 shows X-ray images combining the four 
Chandra observations. Thanks to the long exposure time 
of 190 ks, several fine structures are revealed. The pulsar 
wind nebula is very large, almost filling the radio shell of 
G320.4-01.2. The length of the X-ray nebula along the 
jet-like feature reaches ~ 16', which corresponds to ~ 24 
pc for the distance d = 5.2 kpc. The south edge of the 
nebula is more distant (~ 10') from the pulsar than the 
north edge (~ 6'), while the width of the nebula across 
the east- west direction is about 8', which is a half of the 
length along the axis. 

Interestingly the western side of the PWN appears to 
be sharply cut off at 2' from the pulsar, implying that 
the pulsar wind is confined by an invisible ISM. Moreover 
filament-like absorption features are seen at the western 
edge. Evidently, the column density seems to change at 
these features (Yatsu 2008). A more detailed treatment 
and discussion of these data are deferred to a subsequent 
paper. 

In the close-up image in Figure 1, nested triplet ring 
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features, which are labeled "Outer- Torus" , "Inner- Torus" , 
and "Inner-Ring", are found. The outer features have al- 
ready been reported by Gaensler et al. (2002). The outer 
torus is seen at a separation of ~ 60" from the pulsar. In 
the inner side of the outer-torus, an inner-torus with a 
radius of ~ 30" exists. Most intriguing is an indication 
of an "Inner-Ring" structure ~ 10" apart from the pul- 
sar as seen in the Crab nebula (Weisskopf et al. 2000). 
Nevertheless the existence of the ring feature is still am- 
biguous comparing with the Crab, which might be due to 
the contamination from the bright nebula core. 

3. Data Analysis and Results 

3.1. Imaging analysis 

In order to visualize the compact structures in the vicin- 
ity of the pulsar, we applied an "unsharp mask" , which 
is a kind of the high-pass filter, to the Chandra image. 
First, we generated a low-frequency component image that 
should be removed. To avoid the impact from the bright 
point sources such as the pulsar, we removed noticeable 
point sources using "wavdetect" in CIAO 4.0. The gap 
regions were extrapolated from the surrounding areas us- 
ing "dmfilth" . Then the filled-in image was normalized by 
the exposure map. Finally, the exposure corrected image 
was smoothed using the tool "aconvolve" As a convolu- 
tion kernel, a Gaussian function of a = 3.5" was adopted. 
We also generated a source image by the same method as 
for the low-frequency image but with a smoothing with a 
Gaussian of a = 0.5". Finally, the low-pass filtered com- 
ponent was subtracted from the source image. 

The resultant image is shown in Figure 2. The inner 
ring is clearly distinguished from the pulsar. The radius 
of the ring feature appears to be ~ 10" along the major 
axis, that corresponds to 7.8 x 10 17 cm for the distance of 
5.2 kpc from earth. The eccentricity of e ~ 0.8, for a minor 
axis radius of ~ 6", implies an inclination angle of i ~ 50° 
which is somewhat smaller than that from the appearance 
of the outer torus (Gaensler et al. 2002). On the north 
edge of the inner ring, there arc emitting blobs, which 
forms the north jet. These small structures are seen only 
in the north side of the pulsar, while the south jet appears 
to be diffuse. The south jet across the flow line is about 
20" at the origin, which just corresponds to the radius 
of the ring feature. These structures seems to changes 
temporally (Yatsu 2008). The detailed temporal analysis 
will be reported in a subsequent paper. 

Next we measure the radius of the ring feature based 

Table 1. Summary of the observations of PSR B1509-58. 



ObsID# 


Date 


Detector 


Mode 


Exposure 
(s) 


5534 


2004/12/28 


ACIS-I 


VFaint 


50130 


5535 


2005/02/07 


ACIS-I 


VFaint 


43130 


6116 


2005/04/29 


ACIS-I 


VFaint 


47650 


6117 


2005/10/18 


ACIS-I 


VFaint 


46140 




02.0 15:14:00.0 58.0 56.0 54.0 52.0 13:50.0 48.0 



R.A. (J2000.0) 

Fig. 2. Unsharp masked image of the pulsar vicinity. This 
image was generated from a source image smoothed with a 
Gaussian of a = 0.5", from which a low-pass image smoothed 
with a Gaussian of cr = 3.5" was subtracted. The pulsar's 
location is indicated by the white dashed-line. 

on the radial profile of the surface brightness for a quanti- 
tative discussion. As shown in Figure 3, sixty sub-regions 
were selected from which the X-ray photons were accumu- 
lated. Each of the sub-regions has a width of 1.0". The 
ellipticity of the regions are aligned with the inner-ring. 
In order to exclude the jet components, the north-west 
quadrant and the south-east quadrant respect to the pul- 
sar were masked. Thus the resultant sub-regions exhibit 
bow-tie shapes. 

The radial profile of the X-ray surface brightness for 
0.4^8.0 keV band is shown in Figure 4. The error bars in- 
dicate the statistical uncertainties corresponding the total 
numbers of the X-ray photons detected at the each sub- 
region. Since a sub-region accumulates more than 1000 
photons, the relative uncertainty is smaller than ~ 4 %. 

To constrain the radius of the ring feature, we per- 
formed a model fitting with a multi-component function 
consisting of two power-law functions and a Gaussian 
function. For the fitting we utilized the data points ranged 
in 2" .0 ~ 28" .0. The centroid of the Gaussian that corre- 
sponds to the inner ring was 9".0±0".l at 1 cr confidence 
level. The other parameters are summarized in Table 2. 

3.2. Spectral analysis 

In this section the spatially resolved spectroscopy is pre- 
sented. Since the photon statistics was a little poor, we 
investigated the spectral evolution with a spatial resolu- 
tion of 3" . The examined regions are the north-east side 
of the bow-tie shape in Figure 3. 

For the model fittings, we employed a single compo- 
nent model consisting of an absorbed power-law func- 
tion. However, it is difficult to determine the column 
density for each region individually, because the source 
regions are small. Therefore we fixed the column density 
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Fig. 3. The sub-regions to be analyzed are described by 
white lines. The each sub-region has a width of 1".0 along 
the major axis (North-east^South-wcst). 




Radius [arcsec] 

Fig. 4. The radial profile of the X-ray surface brightness at 
the energy band 0.4-8.0 keV. The solid line shows the model 
function consists of a Gaussian function and two power-law 
functions. The dashed lines show the each component con- 
sisting the model. The dot-dash line shows HRMAs point 
spread function normalized by the pulsar's luminosity at 1" 
from the pulsar. 

at JVh = 1.0 x 10 22 cm -2 . With the above configurations, 
the four data sets observed in the different epochs could 
not be explained by the same parameter set, the energy 
flux varies with time. These results may be caused by 
the temporal changes of the PWN or by the different set- 
tings of instrumentation, such as the roll angle. Thus we 
allowed the normalizations of the four data-sets to vary 
individually. 

The results of model fittings are shown in Figure 
5, indicating the evolution of the surface brightness 
and the photon index as a function of radius. In all 
panels, shapes of the data points correspond to the 
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Fig. 5. The evolutions of the surface brightness (upper 
panel) and the photon index (lower-panel) are shown as func- 
tions of radius. The shape of the data points corresponds 
to the epochs, 2004-Dec-28 (circle), 2005-Feb-07 (triangle), 
2005-Apr-29 (square), and 2005-Oct-18 (star). 

epochs, 2004-Dec-28(circle), 2005-Feb-07(triangle), 2005- 
Apr-29 (square), and 2005-Oct-18(star). The surface 
brightness slightly increases at a radius of ~ 9", corre- 
sponding the inner ring. At R ~ 30", the light curve 
shows a dip which corresponds to the gap between the 
inner torus and the outer-torus. The photon index seems 
to be constant in the inner torus within i?~ 30", while it 
increases in the outer-torus outside the radius of i?~ 30". 

4. Discussion 

4-1. Double tori 

As shown in Figure 4, the PWN around PSR B1509-58 
has a centrally concentrated surface brightness profile, 
which may be caused by its comparatively small angu- 
lar size or the contamination from the bright jet structure 
pointing towards us. For this reason, it is difficult to com- 
pare the surface brightness evolution of PSR B1509— 58 
with that of the Crab nebula in a simple way. Nevertheless 
PSR B1509— 58 possesses an apparently different mor- 
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phology; the nested double ring structure, inner-torus 
(R ~ 30") and outer-torus (30" < R < 60"), lying in the 
equatorial plane relative to the spin axis. 

Spectroscopy, these two tori, calling the inner torus and 
the outer torus, seems to be completely different. The 
inner-torus has a constant photon index of T ~ 1.6 which 
is the same as that just at the pulsar. In contrast, the 
photon index in the outer-torus obviously increases with 
the radius from T = 1.6 to 2.1. The observed interval of 
the photon index variation is Ar ~ 0.5, which implies that 
the observed spectral evolution is caused by synchrotron 
cooling (Kardashev 1962). 

In case of the Crab nebula, the post-shock photon in- 
dex is r = 2.1 and it does not change within the torus. 
According to Kennel and Coroniti (1984a), the evolution 
of the surface brightness in the torus can be explained by 
the compression of the toroidal magnetic field, which en- 
hances the synchrotron energy loss. Thus the bright torus 
of the Crab simply corresponds to the region in which 
most of the X-ray emitting particles radiate their energy 
and burn off there. Consequently the photon index out- 
side of the torus is steeper than T = 2.1, and the surface 
brightness decays rapidly. 

In terms of spectroscopy, the counterpart of the X-ray 
torus in the Crab is the inner torus rather than the outer 
torus of PSR B1509— 58. While in the outer-torus the 
X-ray emitting charged particles have burnt off and thus 
the synchrotron cooling break is passing through the ob- 
served energy band (8.0 keV to 0.4 keV). Moreover, it 
was reported that the outer-torus has a radio counter part 
(Gaensler et al. 2002). These results of the spectroscopy 
imply that the outer-torus could not be simply explained 
by radiative losses as claimed by Gaensler et al. (2002). 

4-2. Inner Ring 

It is essential to know the radius of the pulsar wind 
termination shock to understand the mechanism which 
forms the structure of a PWN. So far, the termination 
shock of the pulsar wind has been observed only in the 
Crab nebula. In this section we discuss the ring feature 
discovered near around PSR B1509— 58 comparing it with 
the inner ring of the Crab nebula. 

As the pulsar wind is believed to be ejected almost at 
the light speed, the free flowing region should not be ra- 
diating as seen in the Crab nebula. In the case of PSR 
B1509— 58 the region within the ring feature is also dark. 
Comparing it with the tori which are strongly affected 
by the Lorcnz boost, the surface brightness of the ring 
feature appears to be uniform, so that we can distinguish 
the south edge of the ring. For the Crab nebula Weisskopf 
et al. 2000 pointed out that the surface brightness of the 
inner ring is more uniform in azimuth, indicating that 
relativistic beaming is less significant in it than in the 
torus. The observed apparent properties can distinguish 
the ring from the tori and then the hypothesis that the 
ring structure represents the termination shock can be 
suggested. By the way, the ring shows an obvious elliptic 
shape, implying that the shock front in an ellipsoidal sur- 
face or a ring viewed with inclination rather than a sphere. 



Regarding the issue, recent numerical simulations predict 
a shock front with an oblate shape for an anisotropic pul- 
sar wind (Komissarov and Lyubarsky 2003; Del Zanna 
et al. 2004). 

At the termination shock, pressure balance between the 
ram pressure of the pulsar wind and the post-shock plasma 
is expected, as seen in the Crab nebula. Based on the X- 
ray image 2, we obtained a radius of the ring structure of 
r ~ 9", which corresponds to 7 x 10 17 cm for the distance of 
5.2 kpc to the pulsar. Assuming that most of the pulsar's 
rotating energy, E = 1.8 X 10 37 ergs s _1 , is transferred to 
the relativistic pulsar wind at the light cylinder, we find 
the ram pressure of the pulsar wind at the ring structure 
to be 



F * 47rr| s 0c y ' 

= 9.6x lO- 11 ^" 1 (—^— 1 (—) 2 dyncm(5) 

where <j) is the fraction of a sphere covered by the wind. 

On the other hand, the inner pressure of the nebula is 
roughly estimated by assuming equipartition conditions 
(Reference). From the spectral fits we obtained an X-ray 
flux of 1.32 x 10 -13 ergs s _1 cm -2 and a photon index 
of r = 1.6 for the region where the ring feature exists 
(R = 10") and we infer an unabsorbed luminosity density 
at 1 keV iikcv = 6.6 x 10 14 ergs s -1 Hz -1 . For sim- 
plification wc assumed the emitting volume to be a par- 
tial spherical shell of radius R = 10", thickness AR = 3" 
and opening angle of 45°, corresponding to V = 2.6 x 10 53 
cm 3 , which may be an overestimate. For a filling fac- 
tor / and a ratio of ion to electron energy /i, the mini- 
mum energy present in the source responsible for the syn- 
chrotron emission Wtotai ~ 2.9(1 + /i) 4 / 7 / 3 / 7 x 10 43 ergs 
and a magnetic field of B ~ 35(1 + /i) 2 / 5 /~ 2 / 7 ^G are in- 
ferred. The corresponding pressure of the relativistic gas 
is then Pneb = Wtotai/3V~ 3.7(1 + /1) 3 / 7 /- 4 / 7 x 10- 11 dyn 
cm~ 2 l . Note that the obtained inner pressure is a lower 
limit because the post-shock flow of a weakly magnetized 
pulsar wind is believed to be in particle dominant condi- 
tions rather than in equipartition (Kennel and Coroniti 
1984a). Taking into account this fact, p* and p nc b are 
comparable, implying that the ram pressure of the pulsar 
wind and the internal pressure of the PWN is balanced 
at the ring feature. It is therefore inferred that the newly 
found ring feature is the termination shock of the pulsar 
wind. 

4-3. Magnetization of the pulsar wind 

Adopting that the discovered ring feature corresponds 
to the wind termination shock, the magnetization param- 
eter a can be estimated. As the shock transforms the 
kinetic energy into thermal energy, a can be rewritten as 
(Shibata et al. 1998) 

1 For a torus geometry with opening angle 15° corresponding to 
/ ~ 0.5, the internal pressure of p ne b ~ 5.6(1 + /j) 3//7 x 10 -11 dyn 
cm -2 can be obtained. 
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[kinetic energy density] 

Bf/An 

[postshock thermal energy density] 
Bf/4n 



(3) 



(4) 



(5) 



where B\ is the toroidal magnetic field just behind the 
shock, and E^q is the equipartition magnetic field after 
the shock. Combining this equation with the equation of 
continuity (r^ s c/3 = r^v n ) and the flux conservation law 
(r TS cBi = r n B Eq ) yields 

T7=- (6) 

Substituting a shock radius tts = 9" and a radius of the 
nebula r n = 30", corresponding to the inner-torus, we ob- 
tain cr ~ 0.01, which is twice larger than that reported by 
Gaensler et al. (2002) although they regarded the inner- 
torus in Figure 1 as a wind termination shock. 

In case of the Crab nebula a magnetization parameter of 
cr = 0.003 based on the expansion velocity was calculated 
by Kennel and Coroniti (1984a), while Mori 2002 also re- 
ported a larger magnetization parameter of ct~ 0.03 based 
on imaging analyses utilizing Chandra. 

4-4- South jet 

The south jet in the pulsar vicinity is remarkably differ- 
ent from that seen in the Crab. In case of the Crab nebula, 
the jet seems to flow out directly from the pulsar, even 
Chandra cannot resolve the origin of the jet (Weisskopf 
et al. 2000; Hester et al. 2002; Mori et al. 2004). In con- 
trast, the jet of PSR B1509-58 extends over ~ 10" at 
the origin. Also, the jet itself is broad and brighter than 
its torus. If the newly discovered ring feature is a wind 
termination shock, the south jet of PSR B1509— 58 may 
flow out from the vicinity of the shock front or from the 
inner-torus, namely the post shock region. 

The appearance of the south jet at the origin seems 
to match the recent theoretical scenario in which the 
hoop stress in the pulsar wind collimatcs the poloidal flow 
(Lyubarsky 2002). This idea was of course introduced for 
the Crab nebula, however, it is difficult to explain the nar- 
row jet because the magnetic pinch works well only down- 
stream of the shock. Downstream, the magnetic hoop 
stress is believed to be able to squeeze the "mildly rel- 
ativistic" flow towards the rotating axis. The observed 
features indicate that this model can be applied to PSR 
B1509-58. 

Del Zanna et al. (2004) reported two dimensional rel- 
ativistic MHD simulations of pulsar jets in various con- 
figurations. Interestingly the observed morphology of the 
PWN around PSR B1509— 58 closely resembles the results 
of the simulations for a = 0.03. They argued that produc- 
ing a jet by magnetic hoop stress requires a comparatively 
large magnetization parameter, a > 0.01, which is coinci- 
dent with the calculated value shown above. It is therefore 
argued that the south jet of PSR B1509— 58 originates in 
the post-shock equatorial flow. 
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Table 2. Summary of the model fitting of the radial profile. 



Parameter Value 

Gaussian 

Normalization 4.76 ±0.44 (xlCT 8 photons s -1 pixel -1 ) 

Mean 9".06±0".14 

Sigma 1".48±0".18 

Power- law (1) 

Normalization 5.36 ±0.65 (xlO -7 photons s _1 pixel -1 ) 

Index -0.69 ±0.04 

Power- law (2) 

Normalization 5.82±1.31 (xlO -6 photons s -1 pixel" 1 ) 

Index -3.48 ±0.37 



The quoted uncertainties are 1 a. 



